Proteogenomic characterization of primary colorectal cancer and metastatic progression identifies proteome-based subtypes and signatures.

Metastatic progression of colorectal adenocarcinoma (CRC) remains poorly understood and poses significant challenges for treatment. To overcome these challenges, we performed multiomics analyses of primary CRC and liver metastases. Genomic alterations, such as structural variants or copy number alterations, were enriched in oncogenes and tumor suppressor genes and increased in metastases. Unsupervised mass spectrometry-based proteomics of 135 primary and 123 metastatic CRCs uncovered distinct proteomic subtypes, three each for primary and metastatic CRCs, respectively. Integrated analyses revealed that hypoxia, stemness, and immune signatures characterize these 6 subtypes. Hypoxic CRC harbors high epithelial-to-mesenchymal transition features and metabolic adaptation. CRC with a stemness signature shows high oncogenic pathway activation and alternative telomere lengthening (ALT) phenotype, especially in metastatic lesions. Tumor microenvironment analysis shows immune evasion via modulation of major histocompatibility complex (MHC) class I/II and antigen processing pathways. This study characterizes both primary and metastatic CRCs and provides a large proteogenomics dataset of metastatic progression.

ATM-deficiency-induced microglial activation promotes neurodegeneration in ataxia-telangiectasia.

While ATM loss of function has long been identified as the genetic cause of ataxia-telangiectasia (A-T), how it leads to selective and progressive degeneration of cerebellar Purkinje and granule neurons remains unclear. ATM expression is enriched in microglia throughout cerebellar development and adulthood. Here, we find evidence of microglial inflammation in the cerebellum of patients with A-T using single-nucleus RNA sequencing. Pseudotime analysis revealed that activation of A-T microglia preceded upregulation of apoptosis-related genes in granule and Purkinje neurons and that microglia exhibited increased neurotoxic cytokine signaling to granule and Purkinje neurons in A-T. To confirm these findings experimentally, we performed transcriptomic profiling of A-T induced pluripotent stem cell (iPSC)-derived microglia, which revealed cell-intrinsic microglial activation of cytokine production and innate immune response pathways compared to controls. Furthermore, A-T microglia co-culture with either control or A-T iPSC-derived neurons was sufficient to induce cytotoxicity. Taken together, these studies reveal that cell-intrinsic microglial activation may promote neurodegeneration in A-T.

Personalized health risk assessment based on single-cell RNA sequencing analysis of a male with 45, X/48, XYYY karyotype.

Numeric sex chromosome abnormalities are commonly associated with an increased cancer risk. Here, we report a 14-year-old boy with a rare mosaic 45, X/48, XYYY karyotype presenting with subtle dysmorphic features and relative height deficiency, requiring growth hormone therapy. As only 12 postnatal cases have been described so far with very limited follow-up data, to assess the proband's long-term prognosis, including cancer risk, we performed high-throughput single-cell RNA sequencing (scRNA-seq) analysis. Although comprehensive cytogenetic analysis showed seemingly near perfect balance between 45, X and 48, XYYY cell populations, scRNA-seq revealed widespread differences in genotype distribution among immune cell fractions, specifically in monocytes, B- and T-cells. These results were confirmed at DNA level by digital-droplet PCR on flow-sorted immune cell types. Furthermore, deregulation of predominantly autosomal genes was observed, including TCL1A overexpression in 45, X B-lymphocytes and other known genes associated with hematological malignancies. Together with the standard hematological results, showing increased fractions of monocytes and CD4+/CD8+T lymphocytes ratio, long-term personalized hemato-oncological surveillance was recommended in the reported patient.

High prevalence of somatic PIK3CA and TP53 pathogenic variants in the normal mammary gland tissue of sporadic breast cancer patients revealed by duplex sequencing.

The mammary gland undergoes hormonally stimulated cycles of proliferation, lactation, and involution. We hypothesized that these factors increase the mutational burden in glandular tissue and may explain high cancer incidence rate in the general population, and recurrent disease. Hence, we investigated the DNA sequence variants in the normal mammary gland, tumor, and peripheral blood from 52 reportedly sporadic breast cancer patients. Targeted resequencing of 542 cancer-associated genes revealed subclonal somatic pathogenic variants of: PIK3CA, TP53, AKT1, MAP3K1, CDH1, RB1, NCOR1, MED12, CBFB, TBX3, and TSHR in the normal mammary gland at considerable allelic frequencies (9 × 10-2- 5.2 × 10-1), indicating clonal expansion. Further evaluation of the frequently damaged PIK3CA and TP53 genes by ultra-sensitive duplex sequencing demonstrated a diversified picture of multiple low-level subclonal (in 10-2-10-4 alleles) hotspot pathogenic variants. Our results raise a question about the oncogenic potential in non-tumorous mammary gland tissue of breast-conserving surgery patients.

Transcriptome Characterization of Matched Primary Breast and Brain Metastatic Tumors to Detect Novel Actionable Targets.

BACKGROUND: Breast cancer brain metastases (BrMs) are defined by complex adaptations to both adjuvant treatment regimens and the brain microenvironment. Consequences of these alterations remain poorly understood, as does their potential for clinical targeting. We utilized genome-wide molecular profiling to identify therapeutic targets acquired in metastatic disease. METHODS: Gene expression profiling of 21 patient-matched primary breast tumors and their associated brain metastases was performed by TrueSeq RNA-sequencing to determine clinically actionable BrM target genes. Identified targets were functionally validated using small molecule inhibitors in a cohort of resected BrM ex vivo explants (n = 4) and in a patient-derived xenograft (PDX) model of BrM. All statistical tests were two-sided. RESULTS: Considerable shifts in breast cancer cell-specific gene expression profiles were observed (1314 genes upregulated in BrM; 1702 genes downregulated in BrM; DESeq; fold change > 1.5, Padj < .05). Subsequent bioinformatic analysis for readily druggable targets revealed recurrent gains in RET expression and human epidermal growth factor receptor 2 (HER2) signaling. Small molecule inhibition of RET and HER2 in ex vivo patient BrM models (n = 4) resulted in statistically significantly reduced proliferation (P < .001 in four of four models). Furthermore, RET and HER2 inhibition in a PDX model of BrM led to a statistically significant antitumor response vs control (n = 4, % tumor growth inhibition [mean difference; SD], anti-RET = 86.3% [1176; 258.3], P < .001; anti-HER2 = 91.2% [1114; 257.9], P < .01). CONCLUSIONS: RNA-seq profiling of longitudinally collected specimens uncovered recurrent gene expression acquisitions in metastatic tumors, distinct from matched primary tumors. Critically, we identify aberrations in key oncogenic pathways and provide functional evidence for their suitability as therapeutic targets. Altogether, this study establishes recurrent, acquired vulnerabilities in BrM that warrant immediate clinical investigation and suggests paired specimen expression profiling as a compelling and underutilized strategy to identify targetable dependencies in advanced cancers.

Immunophenotyping and transcriptional profiling of in vitro cultured human adipose tissue derived stem cells.

Adipose-derived stem cells (ASCs) have become an important research model in regenerative medicine. However, there are controversies regarding the impact of prolonged cell culture on the ASCs phenotype and their differentiation potential. Hence, we studied 10 clinical ASCs replicates from plastic and oncological surgery patients, in six-passage FBS supplemented cultures. We quantified basic mesenchymal cell surface marker transcripts and the encoded proteins after each passage. In parallel, we investigated the differentiation potential of ASCs into chondrocytes, osteocytes and adipocytes. We further determined the effects of FBS supplementation and subsequent deprivation on the whole transcriptome by comprehensive mRNA and miRNA sequencing. Our results show that ASCs maintain differentiation potential and consistent profile of key mesenchymal markers, with apparent expression of distinct isoforms, in long-term cultures. No significant differences were observed between plastic and oncological surgery cohorts. ASCs in FBS supplemented primary cultures are almost committed to mesenchymal lineages as they express key epithelial-mesenchymal transition genes including early mesenchymal markers. Furthermore, combined mRNA/miRNA expression profiling strongly supports a modulatory role for the miR-30 family in the commitment process to mesenchymal lineages. Finally, we propose improvements to existing qPCR based assays that address alternative isoform expression of mesenchymal markers.

Machine learning: a useful radiological adjunct in determination of a newly diagnosed glioma's grade and IDH status.

INTRODUCTION: Machine learning methods have been introduced as a computer aided diagnostic tool, with applications to glioma characterisation on MRI. Such an algorithmic approach may provide a useful adjunct for a rapid and accurate diagnosis of a glioma. The aim of this study is to devise a machine learning algorithm that may be used by radiologists in routine practice to aid diagnosis of both: WHO grade and IDH mutation status in de novo gliomas. METHODS: To evaluate the status quo, we interrogated the accuracy of neuroradiology reports in relation to WHO grade: grade II 96.49% (95% confidence intervals [CI] 0.88, 0.99); III 36.51% (95% CI 0.24, 0.50); IV 72.9% (95% CI 0.67, 0.78). We derived five MRI parameters from the same diagnostic brain scans, in under two minutes per case, and then supplied these data to a random forest algorithm. RESULTS: Machine learning resulted in a high level of accuracy in prediction of tumour grade: grade II/III; area under the receiver operating characteristic curve (AUC) = 98%, sensitivity = 0.82, specificity = 0.94; grade II/IV; AUC = 100%, sensitivity = 1.0, specificity = 1.0; grade III/IV; AUC = 97%, sensitivity = 0.83, specificity = 0.97. Furthermore, machine learning also facilitated the discrimination of IDH status: AUC of 88%, sensitivity = 0.81, specificity = 0.77. CONCLUSIONS: These data demonstrate the ability of machine learning to accurately classify diffuse gliomas by both WHO grade and IDH status from routine MRI alone-without significant image processing, which may facilitate usage as a diagnostic adjunct in clinical practice.

Absence of 1p/19q codeletion in oligodendroglioma-like areas of pilocytic astrocytomas.

The presence of oligodendroglioma-like areas in pilocytic astrocytoma may give rise to pathologic diagnostic uncertainty. This study aims to determine if the oligodendroglioma-like areas present in some pilocytic astrocytomas (PA) possess the signature 1p/19q codeletion that is characteristic of classical oligodendroglioma. Array comparative genomic hybridization was carried out on 12 PA samples, from which oligodendroglioma-like areas were microdissected and used as the template DNA source. 1p/19q codeletions were not found in any of the oligodendroglioma areas in PAs. We conclude that PAs with oligodendroglioma-like areas do not share the same molecular genetics as classic oligodendroglioma.
.

High-Throughput Tabular Data Processor - Platform independent graphical tool for processing large data sets.

High-throughput technologies generate considerable amount of data which often requires bioinformatic expertise to analyze. Here we present High-Throughput Tabular Data Processor (HTDP), a platform independent Java program. HTDP works on any character-delimited column data (e.g. BED, GFF, GTF, PSL, WIG, VCF) from multiple text files and supports merging, filtering and converting of data that is produced in the course of high-throughput experiments. HTDP can also utilize itemized sets of conditions from external files for complex or repetitive filtering/merging tasks. The program is intended to aid global, real-time processing of large data sets using a graphical user interface (GUI). Therefore, no prior expertise in programming, regular expression, or command line usage is required of the user. Additionally, no a priori assumptions are imposed on the internal file composition. We demonstrate the flexibility and potential of HTDP in real-life research tasks including microarray and massively parallel sequencing, i.e. identification of disease predisposing variants in the next generation sequencing data as well as comprehensive concurrent analysis of microarray and sequencing results. We also show the utility of HTDP in technical tasks including data merge, reduction and filtering with external criteria files. HTDP was developed to address functionality that is missing or rudimentary in other GUI software for processing character-delimited column data from high-throughput technologies. Flexibility, in terms of input file handling, provides long term potential functionality in high-throughput analysis pipelines, as the program is not limited by the currently existing applications and data formats. HTDP is available as the Open Source software (https://github.com/pmadanecki/htdp).

The DNA copy number landscape of a collision tumor.

Intracranial collision tumors are composed of two histologically distinct but merging components, and are rare. Their genetic profile has rarely been described. Comparative genome hybridization of a combined meningioma and oligodendroglioma demonstrated deletion of chromosome 22q and of 19q in both tumors. Somatic deletion of chromosome 22q and 19q is associated with development of an intracranial collision tumor.
.

Concurrent DNA Copy-Number Alterations and Mutations in Genes Related to Maintenance of Genome Stability in Uninvolved Mammary Glandular Tissue from Breast Cancer Patients.

Somatic mosaicism for DNA copy-number alterations (SMC-CNAs) is defined as gain or loss of chromosomal segments in somatic cells within a single organism. As cells harboring SMC-CNAs can undergo clonal expansion, it has been proposed that SMC-CNAs may contribute to the predisposition of these cells to genetic disease including cancer. Herein, the gross genomic alterations (>500 kbp) were characterized in uninvolved mammary glandular tissue from 59 breast cancer patients and matched samples of primary tumors and lymph node metastases. Array-based comparative genomic hybridization showed 10% (6/59) of patients harbored one to 359 large SMC-CNAs (mean: 1,328 kbp; median: 961 kbp) in a substantial portion of glandular tissue cells, distal from the primary tumor site. SMC-CNAs were partially recurrent in tumors, albeit with considerable contribution of stochastic SMC-CNAs indicating genomic destabilization. Targeted resequencing of 301 known predisposition and somatic driver loci revealed mutations and rare variants in genes related to maintenance of genomic integrity: BRCA1 (p.Gln1756Profs*74, p.Arg504Cys), BRCA2 (p.Asn3124Ile), NCOR1 (p.Pro1570Glnfs*45), PALB2 (p.Ser500Pro), and TP53 (p.Arg306*). Co-occurrence of gross SMC-CNAs along with point mutations or rare variants in genes responsible for safeguarding genomic integrity highlights the temporal and spatial neoplastic potential of uninvolved glandular tissue in breast cancer patients.

Modulation of chemotherapeutic drug resistance in neuroblastoma SK-N-AS cells by the neural apoptosis inhibitory protein and miR-520f.

The acquisition of multidrug resistance is a major impediment to the successful treatment of neuroblastoma, a clinically heterogeneous cancer accounting for ∼15% of all pediatric cancer deaths. The MYCN transcription factor, whose gene is amplified in ∼30% of high-risk neuroblastoma cases, influences drug resistance by regulating a cadre of genes, including those involved with drug efflux, however, other high-risk subtypes of neuroblastoma lacking MYCN amplification, such as those with chromosome 11q deletions, also acquire multidrug resistance. To elucidate additional mechanisms involved with drug resistance in non-MYCN amplified tumour cells, an SK-N-AS subline (SK-N-AsCis24) that is significantly resistant to cisplatin and cross resistant to etoposide was developed through a pulse-selection process. High resolution aCGH analysis of SK-N-AsCis24 revealed a focal gain on chromosome 5 containing the coding sequence for the neural apoptosis inhibitory protein (NAIP). Significant overexpression of NAIP mRNA and protein was documented, while experimental modulation of NAIP levels in both SK-N-AsCis24 and in parental SK-N-AS cells confirmed that NAIP was responsible for the drug resistant phenotype by apoptosis inhibition. Furthermore, a decrease in the NAIP targeting microRNA, miR-520f, was also demonstrated to be partially responsible for increased NAIP levels in SK-N-AsCis24. Interestingly, miR-520f levels were determined to be significantly lower in postchemotherapy treatment tumours relative to matched prechemotherapy samples, consistent with a role for this miRNA in the acquisition of drug resistance in vivo, potentially through decreased NAIP targeting. Our findings provide biological novel insight into neuroblastoma drug-resistance and have implications for future therapeutic research.

miR-199a-5p silencing regulates the unfolded protein response in chronic obstructive pulmonary disease and α1-antitrypsin deficiency.

RATIONALE: Retention of abnormal α1-antitrypsin (AAT) activates the unfolded protein response in AAT-deficient monocytes. The regulatory role of microRNAs (miRNAs) in unfolded protein responses and chronic obstructive pulmonary disease pathogenesis has not been investigated. OBJECTIVES: To investigate miRNA expression and function in MM and ZZ monocytes and identify miRNA(s) regulating the unfolded protein response. METHODS: Peripheral blood monocytes were isolated from asymptomatic and symptomatic MM and ZZ individuals for miRNA expression profiling and pyrosequencing analysis. miRNA/gene and protein expression was measured with quantitative polymerase chain reaction and Western blotting. Overexpression and inhibition studies were performed with pre-miR or anti-miR, respectively. Luciferase reporter genes were used to elucidate direct miRNA-target interactions. Inflammatory cytokines were detected using the Meso Scale Discovery Plex assays. MEASUREMENTS AND MAIN RESULTS: Forty-three miRNAs were differentially expressed, with miR-199a-5p most highly up-regulated in asymptomatic ZZ versus MM monocytes. miR-199a-2 promoter hypermethylation inhibits miR-199a-5p expression and was increased in symptomatic MM and ZZ monocytes compared with asymptomatic counterparts. GRP78, activating transcription factor 6, p50, and p65 were increased in symptomatic versus asymptomatic ZZ monocytes. Reciprocal down- or up-regulation of these markers was observed after miRNA modulation. Direct miR-199a-5p targeting of activating transcription factor 6, p50, and p65 by miR-199a-5p was demonstrated using luciferase reporter systems. Overexpression of miR-199a-5p also decreased other arms of the UPR and expression of cytokines that are not putative targets. CONCLUSIONS: miR-199a-5p is a key regulator of the unfolded protein response in AAT-deficient monocytes, and epigenetic silencing of its expression regulates this process in chronic obstructive pulmonary disease.

Regeneration of comparative genomic hybridization oligonucleotide microarrays with dimethylurea.

Reuse of materials in DNA hybridization-based methods has been known since the advent of Southern membranes. Array-based comparative genomic hybridization is essentially Southern hybridization with multiple probes immobilized on a solid surface. We show that comparative genomic hybridization microarrays fabricated with maskless array synthesizer technology can be used up to four times with the application of 1,3-dimethylurea as an array-stripping agent. We reproducibly detected chromosomal aberrations (0.6-22.4Mb in size) in four hybridization rounds using regenerated microarray slides. We also demonstrated that regenerated arrays can detect smaller alterations (16-200kbp), such as common copy number variants, as well as complex aberration profiles in tumor DNA.

Differential DNA methylation patterns define status epilepticus and epileptic tolerance.

Prolonged seizures (status epilepticus) produce pathophysiological changes in the hippocampus that are associated with large-scale, wide-ranging changes in gene expression. Epileptic tolerance is an endogenous program of cell protection that can be activated in the brain by previous exposure to a non-harmful seizure episode before status epilepticus. A major transcriptional feature of tolerance is gene downregulation. Here, through methylation analysis of 34,143 discrete loci representing all annotated CpG islands and promoter regions in the mouse genome, we report the genome-wide DNA methylation changes in the hippocampus after status epilepticus and epileptic tolerance in adult mice. A total of 321 genes showed altered DNA methylation after status epilepticus alone or status epilepticus that followed seizure preconditioning, with >90% of the promoters of these genes undergoing hypomethylation. These profiles included genes not previously associated with epilepsy, such as the polycomb gene Phc2. Differential methylation events generally occurred throughout the genome without bias for a particular chromosomal region, with the exception of a small region of chromosome 4, which was significantly overrepresented with genes hypomethylated after status epilepticus. Surprisingly, only few genes displayed differential hypermethylation in epileptic tolerance. Nevertheless, gene ontology analysis emphasized the majority of differential methylation events between the groups occurred in genes associated with nuclear functions, such as DNA binding and transcriptional regulation. The present study reports select, genome-wide DNA methylation changes after status epilepticus and in epileptic tolerance, which may contribute to regulating the gene expression environment of the seizure-damaged hippocampus.

Identification of a novel recurrent 1q42.2-1qter deletion in high risk MYCN single copy 11q deleted neuroblastomas.

Neuroblastoma is an aggressive embryonal tumor that accounts for ∼15% of childhood cancer deaths. Hitherto, despite the availability of comprehensive genomic data on DNA copy number changes in neuroblastoma, relatively little is known about the genes driving neuroblastoma tumorigenesis. In this study, high resolution array comparative genome hybridization (CGH) was performed on 188 primary neuroblastoma tumors and 33 neuroblastoma cell lines to search for previously undetected recurrent DNA copy number gains and losses. A new recurrent distal chromosome 1q deletion (del(1)(q42.2qter)) was detected in seven cases. Further analysis of available array CGH datasets revealed 13 additional similar distal 1q deletions. The majority of all detected 1q deletions was found in high risk 11q deleted tumors without MYCN amplification (Fisher exact test p = 5.61 × 10(-5) ). Using ultra-high resolution (∼115 bp resolution) custom arrays covering the breakpoints on 1q for 11 samples, clustering of nine breakpoints was observed within a 12.5-kb region, of which eight were found in a 7-kb copy number variable region, whereas the remaining two breakpoints were colocated 1.4-Mb proximal. The commonly deleted region contains one miRNA (hsa-mir-1537), four transcribed ultra conserved region elements (uc.43-uc.46) and 130 protein coding genes including at least two bona fide tumor suppressor genes, EGLN1 (or PHD2) and FH. This finding further contributes to the delineation of the genomic profile of aggressive neuroblastoma, offers perspectives for the identification of genes contributing to the disease phenotype and may be relevant in the light of assessment of response to new molecular treatments.

Co-localization of the oncogenic transcription factor MYCN and the DNA methyl binding protein MeCP2 at genomic sites in neuroblastoma.

BACKGROUND: MYCN is a transcription factor that is expressed during the development of the neural crest and its dysregulation plays a major role in the pathogenesis of pediatric cancers such as neuroblastoma, medulloblastoma and rhabdomyosarcoma. MeCP2 is a CpG methyl binding protein which has been associated with a number of cancers and developmental disorders, particularly Rett syndrome. METHODS AND FINDINGS: Using an integrative global genomics approach involving chromatin immunoprecipitation applied to microarrays, we have determined that MYCN and MeCP2 co-localize to gene promoter regions, as well as inter/intragenic sites, within the neuroblastoma genome (MYCN amplified Kelly cells) at high frequency (70.2% of MYCN sites were also positive for MeCP2). Intriguingly, the frequency of co-localization was significantly less at promoter regions exhibiting substantial hypermethylation (8.7%), as determined by methylated DNA immunoprecipitation (MeDIP) applied to the same microarrays. Co-immunoprecipitation of MYCN using an anti-MeCP2 antibody indicated that a MYCN/MeCP2 interaction occurs at protein level. mRNA expression profiling revealed that the median expression of genes with promoters bound by MYCN was significantly higher than for genes bound by MeCP2, and that genes bound by both proteins had intermediate expression. Pathway analysis was carried out for genes bound by MYCN, MeCP2 or MYCN/MeCP2, revealing higher order functions. CONCLUSIONS: Our results indicate that MYCN and MeCP2 protein interact and co-localize to similar genomic sites at very high frequency, and that the patterns of binding of these proteins can be associated with significant differences in transcriptional activity. Although it is not yet known if this interaction contributes to neuroblastoma disease pathogenesis, it is intriguing that the interaction occurs at the promoter regions of several genes important for the development of neuroblastoma, including ALK, AURKA and BDNF.

MicroRNA and DNA methylation alterations mediating retinoic acid induced neuroblastoma cell differentiation.

Many neuroblastoma cell lines can be induced to differentiate into a mature neuronal cell type with retinoic acid and other compounds, providing an important model system for elucidating signalling pathways involved in this highly complex process. Recently, it has become apparent that miRNAs, which act as regulators of gene expression at a post-transcriptional level, are differentially expressed in differentiating cells and play important roles governing many aspects of this process. This includes the down-regulation of DNA methyltransferases that cause the de-methylation and transcriptional activation of numerous protein coding gene sequences. The purpose of this article is to review involvement of miRNAs and DNA methylation alterations in the process of neuroblastoma cell differentiation. A thorough understanding of miRNA and genetic pathways regulating neuroblastoma cell differentiation potentially could lead to targeted therapies for this disease.